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We propose a new method for investigating the global properties of the retarded Green's function 
Gr{x' ,x) for fields propagating on an arbitrary globally hyperbolic spacetime. Our method combines 
the Hadamard form for Gr (this form is only valid within a normal neighbourhood of a;) together 
with Kirchhoff 's integral representation for the field in order to calculate Gr outside the maximal 
normal neighbourhood of a;. As an example, we apply this method to the case of a scalar field on a 
black hole toy-model spacetime, the Plebahski-Hacyan spacetime, M2 x The method allows us 
to determine in an exact manner that the singularity structure of the 'direct' term in the Hadamard 
form for Gr{x',x) changes from a form 5{a) to 'l/yra' after the null geodesic joining x and x' has 
crossed a caustic point, where a is the world function. Furthermore, there is a change of form 
from a 9{—a) to a ln|(T|/7r' in the 'tail' term, which has not been explicitly noted before in the 
literature. We complement the results from the Kirchhoff integral method with an analysis for 
large-^ of the Green function modes. This analysis allows us to determine the singularity structure 
after null geodesies have crossed an arbitrary number of caustics, although it raises a causality issue 
which the Kirchhoff integral method resolves. Because of the similarity in the caustic structure of 
the spacetimes, we expect our main results for wave propagation to also be valid on Schwarzschild 
spacetime. 



I. INTRODUCTION 



Understanding the propagation of waves on a curved spacetime is of interest on its own right, and also has many 
important applications in various physical settings. Consider, for example, the inspiral of a black hole binary in the 
extreme mass ratio limit. This system can be modelled by the motion of the smaller black hole in the gravitational field 
of the larger black hole (typically Schwarzschild or Kerr black hole) under the action of a self-force P] induced by its 
own gravitational field. At the heart of the analysis of such systems are the gravitational wave perturbations generated 
by the small black hole in the curved spacetime of the large black hole. In particular, the self-force is governed by 
an equation which involves the integral of the gradient of the retarded Green's function of the corresponding wave 
equation. The integral is taken over the entire past of the world line of the particle. While the local structure of 
the Green's function is well known in terms of the so-called Hadamard form pj, a lot less is known about the global 
structure of the Green's function. Clearly, however, the structure - in particular, the singularity structure - of the 
Green's function for spacetime points arbitrarily separated is crucial for the understanding of wave propagation and 
it can be of practical use for the calculation of the self-force. This is the basic aim of the present paper: to introduce 
a method that allows the calculation of the retarded Green's function at arbitrarily separated spacetime points. We 
apply this method to determine new results regarding the singularity structure of the retarded Green's function in a 
black hole toy-model spacetime. 

It is known that the singularities of the Green's function occur whenever the two spacetime points are connected 
by a null geodesic (see e.g. Oil])- Within a normal neighbourhood of the base point x (i.e. a region N{x) containing 
X with the property that every x' £ Af{x) is connected to a; by a unique geodesic which lies in Af{x)), the Hadamard 
form dictates that the singularity of the retarded Green function is of the type S{a), where a = a(x,x') is Synge's 
world- function, i.e., one half of the geodesic distance along the unique geodesic connecting the points x and x' . In 
principle, outside a normal neighbourhood the biscalar a{x, x') is not well-defined and the Hamadard form is not valid; 
therefore, one must take an alternative approach to the calculation of the Green function. In 5 it was shown by 
using a quasinormal mode series that the scalar Green's function Gji{x,x') presents a four-fold singularity structure: 
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(5(cr), P.V. (I/ttct) , —6{a), — P.V. (I/ttct) , S{a), . . . , where the change of character in the singularity occurs every time 
that the null geodesic connecting the two points passes through a caustic (i.e., a spacetime point where neighbouring 
null geodesies are focused). This was shown specifically for the case of a static region of Nariai spacetime {682 x S^) but 
it was anticipated to be also valid in most situations in spherically symmetric spacetimes. (We anticipate exceptional 
points in every spherically symmetric spacetime - see Eq.(131) - and indeed there are exceptional spacetimes, for 
example Bertotti-Robinson spacetime.) Such a four-fold structure in General Relativity had been previously noted by 
Ori [B]. It was subsequently shown in 7], again by using a quasinormal mode series, that this four- fold structure is also 
present in Schwarzschild spacetime - although, in this case, the singularity times given by a (truncated) quasinormal 
mode series only approximate the singularity times given by a null geodesic connecting the two spacetime points. 

Furthermore, very recently it has been shown in [8] that it is possible to probe the global singularity structure of the 
Green's function on a generic spacetime by exploiting the Penrose limit [5] and results relating singularity structures 
and caustics in plane wave spacetimes. In particular, the authors find that the four-fold singularity structure applies 
in a wide class of spacetimes including Schwarzschild (and most likely, Kerr). 

A causality issue is raised by the presence of the singular term P.V. (I/ttct) in its naive form indicated above: it 
raises a question regarding the causal character of the Green's fimction, since P.V. (I/ttct) is nonzero even when the 
spacetime points are spacelike separated. While the quasinormal mode method is of practical application (albeit 
to special choices of the mass and, in non-Ricci-fiat spacetimes, of the coupling constant), it does not answer this 
question and it does not generally offer exact control of the terms neglected. The same holds for the Penrose limit 
approach of [5]. This issue is resolved in our approach, which also has the advantage of deriving a corresponding 
four-fold structure for the so-called tail term of the retarded Green's function. We note in particular that this term 
contributes to the singular part of the Green's function after the formation of (an odd number of) caustics. 

The normal neighbourhood, and hence also the Hadamard form for Gn, generally breaks down at caustic points 
of the background spacetime. In a spherically symmetric spacetime, the spray of null geodesies emitted from an 
event x can include an S^-envelope that reconverges at a point x' - a caustic - whose angular separation from x is 
7 = TT. It seems to be precisely the focusing of this S^-envelope of null geodesies that yields the four-fold singularity 
structure in Gji{x, x') for points with 7 ^ 0, tt. (For example, in the case of the Einstein Static Universe, M x S'^, the 
singularity structure is two-fold, and of course the singularity structure is uniform in fiat spacetime.) This is the case 
in Schwarzschild spacetime, but it is also the case in a spacetime M — AI2 x S^, the direct product of a 2-dimensional 
spacetime M2 of constant curvature and the unit 2-sphere (both equipped with their standard metrics). Such 
spacetimes have found application in different contexts in general relativity (see section 7.2 of [10^)- 

In this paper we develop a new method for the formal construction of the retarded Green's function that applies 
beyond the normal neighbourhood. Consider a field 4> propagating on a globally hyperbolic spacetime, whose evolution 
is governed by a homogeneous wave equation. Kirchhoff's formula [11 gives a representation of (/> at a certain point x 
as a convolution integral of G^, the retarded Green's function of the wave equation, and its derivative, with initial 
data for (f> on a spacelike hypersurface in the past of x. In the Kirchhoff integral, Gr acts as a propagator for the 
initial data for (f). By noting that Gii{x,x') satisfies the homogeneous wave equation if x' lies strictly to the future 
of X, we can apply Kirchhoff's formula to the retarded Green's function itself instead of the field. In our method, 
which we refer to as the Kirchhoff integral method, we apply Kirchhoff's formula with acting both as propagator 
and data. We choose time separations so that the Hamadard form is valid for both versions of Gr appearing in the 
Kirchhoff (propagator and data): this puts an upper bound on the time step we can take. However, through repeated 
application of Kirchhoff's formula in this manner we can, in theory, calculate Gn for points arbitrarily separated with 
the sole knowledge of the Hadamard form for Gr. 

We use this method to study the singularity structure of the Green's function as the propagating field passes 
through caustic points of the background spacetime. We note that in [11] a method was employed which has certain 
similarities with our Kirchhoff integral method, although they are fundamentally different methods. The Kirchhoff 
integral method also essentially underpins the standard derivation of the reciprocity relationship Gr{x, x') = Ga{x' , x) 
between the retarded and advanced Green's functions dHH], and to determine properties of the Detweiler- Whiting 
singular Green's function: see Section 14.5 of [1]. Our method is in principle valid for any globally-hyperbolic 
spacetime. However, in practise one must have a knowledge of the Hadamard form for Gr. 

In this paper, we apply the Kirchhoff integral method to the case of a scalar field propagating in a black hole 
toy-model spacetime M — M2 x S^, where M2 is 2-dimensional Minkowski spacetime. Following [15], we will refer 
to this spacetime M as Plebahski-Hacyan (PH) spacetime. We investigate properties of wave propagation on PH 
in the case of a scalar field created by a scalar charge. However, many of the properties of wave propagation are 
independent of the spin of the propagating field and we believe that the approach presented here would also be valid 
in the higher-spin cases. 

Given that the calculations are rather involved, we present here the main result that we obtain when applying the 
Kirchhoff integral method to PH. For field points x' in a normal neighbourhood of the base point x, the Hadamard 
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form for the retarded Green's function applies. This is given by [TJ [31 [T3] 

Gr{x,x') = [U{x,x')S{a{x,x')) + V{x,x')H{-<7{x,x'))]H+ix,x'), (1) 
where S, H are the usual Dirac delta and Heaviside distributions, 



H+{x,x) 



1 if a; lies to the future of : 
otherwise, 



and U,V are biscalars determined as follows. It can be shown that U{x,x') = A^^^(x, x'), where A{x,x') is the 
van Vleck-Morette determinant (which again, is defined only when x' lies in a normal neighbourhood of x). The 
bitensor V{x, x') can be written as an asymptotic series, V{x, x') = X^^o ^nix, a:')cr", where x') are coefficients 
satisfying certain recurrence relations [HIH]. We refer to the terms US{a) and VH{—<t) in the Hadamard form ([T]) 
as, respectively, the 'direct' and the 'tail' terms. 

The world function of PH spacetime can be written as a — —\if' + \l'^ where 77 is the geodesic distance on M2 and 
7 is the geodesic distance on S^. Then for points x' in a normal neighbourhood of x, the Hadamard form applies with 

1/2 00 

, V{x,x') = y^yk{l)<y''. 



U{x,x') = U{j) = 



7 



sm7 



k=0 



where Vk are smooth functions of 7 satisfying certain recurrence relations. Notice then that we can write V{x,x') — 
T(77,7). The normal neighbourhood condition breaks down at 77 = vr, corresponding to the formation of the first 
caustic. The second caustic forms at r/ = 2tt. Our main result from the application of the Kirchhoff integral method 
to PH is that, for tt < 77 < 2tt, 

Gn{x,x')^U{27r~j)Pv( - ^o(27r - 7)^^^^^^^ + (2) 



7r(Ti(x, x') 



where 



a^ = -lv' + li2n-if, (3) 
and the 0(1) term corresponds to a function that remains finite throughout the region tt < rj < 2tt. We note that 

vo{2tt~j) = T(r/,27r-7)|„^=o- 

That is, we show that, after the wavefront has crossed one caustic, the 'direct' term in the Green's function changes 
its singularity character from 6{(t) to P.V. (I/ttct), as expected, and that the 'tail' term gives rise to a singular term 
In I cr I in the same region. This provides the complete description of the singular part of the retarded Green's function 
after the wavefront has crossed one caustic point. To the best of our knowledge, the diverging term ln|cri| has not 
been noted before explicitly in any spacetime. We expect a similar singularity to occur in Schwarzschild spacetime. 

We note that the four-fold singularity structure degenerates to a two-fold singularity structure at the caustic point 
itself, which entails spacetime points with angular separation 7 = tt. This is demonstrated explicitly below: see 



Section IV-F and especially Eq.( 131 ). The dominant part of the singularity displays the form 6{a) — —S{cr) at 7 = tt. 

Our method answers cleanly the two questions raised above, the one about an extension of the world function a 
outside the normal neighbourhood (here given by CTi in (|3])) and the one about 'acausality' due to P.V. (l/Tra) (the 
result ([2]) here is only valid for tt < 77 < 27r), by treating directly the spacetime in question and obtaining exact results. 
Furthermore, we clearly identify the coefficients of the PV and In singularities, suggestively expressed in terms of the 
bitensors U and 1^0 defined in a normal neighbourhood of x. 

We complement the new Kirchhoff integral method with a large-£ asymptotic analysis of the Green's function modes 
in PH spacetime, where £ is the multipole number. This analysis is very similar in nature to the quasinormal mode 
method used in [SJ [7] (although here it is not necessary to perform a Fourier transform and so the modes correspond to 
a multipole decomposition only). As such, and unlike the Kirchhoff integral method, the large--^ asymptotic analysis 
uncovers the singularity structure for an arbitrary number of caustic crossings, but it does not resolve the causality issue 
regarding the term P.V. (I/ttct) mentioned above. The large-^ analysis that we carry out here yields not only the four- 
fold structure for the 'direct' term {±S{a) — >■ =fP-V. (I/ttct)) but also that for the 'tail' term {±6{—a) — >■ =Flii|cr| /tt), 
in both cases for any number of caustic crossings. FigjS] shows the four-fold structure as obtained from an £-mode 
series. A shortcoming of the large-£ asymptotic analysis is that we must restrict to a particular value of + S^R, 
where m is the mass, ^ is the coupling term and R is the Ricci scalar. 
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The layout of the paper is as follows. In Secffl we describe the Kirchhoff integral method as can be applied to a 
general, globally-hyperbolic spacetime. In Sec III we consider briefly the world function and the Hadamard form for 
Gft, in spacetimes of the type Mi x M2 and, in particular, in PH spacetime. In Sec |IV| we apply the Kirchhoff integral 
method specifically to the PH spacetime and provide the lengthy and detailed calculations leading to the result 
This long section is almost entirely technical in nature, and stands independent of the remainder of the paper except 
insofar as it provides the proof of our main result. The large-^ calculation is given in Sec|V]and we conclude with 
some final remarks in SeclVIl 



II. KIRCHHOFF METHOD FOR THE RETARDED GREEN'S FUNCTION 



We consider the wave equation for a scalar field '^(x), 

(□ - - = 0, (4) 

where □ is the d'Alembertian operator, m is the scalar mass, R is the Ricci scalar and ^ is a coupling constant. The 
retarded Green's function for this equation satisfies 

(□ - - ^R)Gr{x, x') = -AnSiix, x'), (5) 

along with the boundary conditions Gret{x" ,x) = if x" ^ J^{x). The quantity 84 is the 4-dimensional Dirac 
distribution. A useful property of Gr is that it allows us to explicitly write down the value of a field ^' at a point x 
to the future of an initial data hypersurface S' [1] [2] : 

^{x) ^ JjGn{x,x')W''' ^{x') ~ ^{x')W''' Gnix,x'))d^'^,, (6) 

where dT,'^, = —ria'dV' is the surface element on E', the future directed normal one-form and dV' is the invariant 
volume element on E'. Thus the Kirchhoff integral Q involves the convolution of Gr and its derivative with the 
initial data for on E'. 

Consider now a globally hyperbolic spacetime {M^gap) with a global time function T. Then the hypersurfaces 

Et = {a; e Af : T{x) = t] 

are Cauchy surfaces for the spacetime. At this point, we implement the following notational change. We take x" to 
be the base point for our calculations, and consider points on hypersurfaces E^/ and Ej that lie to the future of the 
hypersurface Ej//. So with 

T{x) = t > T[x') ^t' > T(x") = t\ 

we see that the retarded Green's function Gfi{x,x") satisfies the homogeneous wave equation in the x variables, and 
so we can apply the Kirchhoff integral ^ to write down 

Gr{x, x") = ^ / {Gr{x, x')dt'GR(x', x") - dt'GRix, x')Gr{x', x"))dV'. (7) 

In this integral, Gr{x,x') acts as the propagator for the field Gr{x' , x"), carrying it from E^ to E^. We note that 
([7]) is well-defined as a distribution, being the convolution of distributions with compact support: see sections 2.5 and 
2.9 of 13]. 

The formula ([7| is of particular interest when the time intervals t' — t" and t — t' are such that both Gr{x',x") 
and Gr{x,x') may be written in Hadamard form, but the field point x does not lie in a normal neighbourhood of 
the base point x" . Hence ([t]) (formally) provides a closed form for the retarded Green's function beyond the normal 
neighbourhood. By applying ([7| iteratively in this way, we have a mechanism for calculating Gr globally. 



HI. PRODUCT SPACETIMES Mi x M2 



General case 



Let us consider a product manifold M = Mi x M2, i.e., ds'^ = dsf -I- ds^ where dsf = g^sdx^dx^ and = 
gabdx'^dx^ . Greek letters a, /3 run over all the coordinates in M, capital Latin letters A, B run over the coordinates in 
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Ml and small Latin letters a, b over the coordinates in M2. Within this subsection only, subindices 1 and 2 will indicate 
quantities corresponding to the manifolds Mi and M2, respectively. Then the defining property VqCtV^it — 2a implies 
that 



cr = (Ti + a-2, 
VaCTiV^cti = 2(71, 

V„(72 W2 = 2(72 



(8) 



We now assume that Mi is a 2-dimensional spacetime, M2 is a 2-dimensional Riemannian manifold and so M is a 
4-dimensional spacetime. Defining A and fj, via ai — — A^/2 and (72 = M^/2 we then have that 



These immediately imply that 



VaAV^A = -1, □iCTi = l-AniA 
S/alJiytJ- = +1, 02^2 = 1 + ^J!^2^^ 



OF dF 
VW(A,M) = A-+,- 



(9) 



(10) 



(9A2 d\ 



for a sufficiently smooth function F of A and /i. From the equation satisfied by the van Vleck determinant A(x,.t'), 
(Act") = 4A, we find 



A = Ai-A2, 
a^VA^i + AiDio-i = 2Ai 

VaA2 + A2n20-2 = 2A2 



(11) 



As shown in jjj, the restriction of the Hadamard biscalar V{XtX') to the null cone V :— V\cr=vi — t'o|cr=o satisfies 
a certain transport equation along the null geodesies ruling the null cone. The bitensor uq{x^x') is defined by 
V ~ Yl^=o '^ncr", as given in the introduction. We can take A~^/^i^o ~ xo{X) + 2/o(a*) to obtain 



- /i(A) + l{m^+m = - /2(m) - iK = const. 



(12) 



2/i(A)-A- 



- A-1/2 



''^i 'DiA 1 



rfA2 



dA 



2/2 (a^) = A^'^' 



-□2A* 



dA 



1/2 



where fi and /2 are some functions of A and ^, respectively. It then follows that 



= + / dA /i(A) + / /2(m) + ^ + 



A 



(13) 



where c^o and Cyo are constants of integration to be determined by the condition that [vq] — [V] = (1 — 3(to^ + 
^i?(a;))/6, where the brackets [ ] around a certain bitensor A{x,x') indicate the coincidence limit of that bitensor (i.e., 
limj;_5.j;/ A{x,x')). This condition follows from the transport equation that V satisfies. 



B. PH: 



X §^ 



The line element of PH spacetime can be written as 

= -dT^ + dY^ + dfl^ , 

where (T, Y) are global inertial coordinates on M2 and dfl^ is the standard line element on : 

dn^ = de^ + sin^ ed$2 

in standard polar and azimuthal coordinates (0, $). The Ricci tensor Rap and Ricci scalar satisfy 

Rcpdx°'dx^ = dn^, R = 2 



(14) 
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and the only non-vanishing Weyl curvature component can be written in Newman-Penrose form as 

1 

'6' 



1p2 



The d'Alembertian is given by 
We find that 



□S2 = 



cot 7 



cr{x, x') 



1 



1 



d_ _ 1 d 
97 sin 7 (?7 

1 2 
o7i 



sm7 



57 



where t' = T(a;'), 6* = 6(0;) etc and 

cos 71 = cos 6* cos 9' + sin 6* sin 9' cos( 
We write the world function on M2 as 



(j{x, x') 



(15) 

(16) 
(17) 
(18) 



We note that when x' lies i n a normal neighbourhood of x, a as given here is the world function of PH spacetime, 
but we will also understand (16) to serve as the global definition of a two point function on M2 x S^. 

In this spacetime, it is straightforward to show that the §~-envelope of null geodesies from any point x" with 
T{x") = t" reconverges at a point x' with T{x') — t" + n. If T{x') < T{x") + tt, then x" and x' are connected by at 
most one causal geodesic. Hence, in line with our comment at the end of Sec|llj the Kirchhoff formula ([T]) is of most 
interest when 



and when 



< i' - t" < TT and < i - t' < TT 



Ti <t-t" < 27r. 



Fig[T] offers a visualization of null geodesies propagating on M2 x §2 on two different T = const, slices. 
We consider now the Hadamard form on PH spacetime. We can apply (11) to obtain 



U{x,x')^U{j) 



sm7 



1/2 



(19) 
(20) 

(21) 



In a general spacetime, the coefficient V{x,x') of the tail term in ([TJ can be written in the form V = J2'^=o ^nO""- 
This is an asymptotic series for V which converges uniformly in a normal neighbourhood of a; in a general class 
of spacetimes [21 [T3]. The coefficients are obtained by solving certain recurrence relations, which correspond to a 
sequence of transport equations pi I14j. In PH spacetime, we find that = i^ki'j) so that 



k=0 



and the transport equations become ordinary differential equations. Applying Eq.(13), we obtain 



V = V\,^o = Ml) - o ^^(7) 1 - 4m^ - + 



r 



cot 7 

7 



Writing h'k — U{"i)vk{l), the recurrence relations become 



(^7 



1 7' 



fc-1/2 



2k \/sm^ I (i7 



sm7 



d 



7 



^sin7 



1/2 



(m^ + ^i?) sin 7 



sm7 
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From Eqs.(23) and (24) we find 



272 - 3csc2(7) [672 + 27sin(27) + 5cos(27) - 5] 



(22) 



(23) 



(24) 



(25) 



25674 

for the specific case + £,R =1/4 and where we choose the constant of integration so that i>i is regular at 7 = 

Unfortunately, we are not able to go any further and find 1^2 from ( 24 ) . 

In Figj3] we compare V{x, x') and vq + vicr with the exact Gji{x,x') as obtained using an ^-mode sum (see Sec V| 
In the next section we evaluate ^ in the PH spacetime and in Sec VI we discuss the consequences of the result. 



FIG. 1: (Colour online.) Visualization of PH spacetime M2 x at given spacelike slices T = To = constant: To € (0, tt) in the 
top figure (referred to as (a) below) and To € (tt, 27r) in the bottom figure (b). The coordinate along the axis of the cylinder 
corresponds to Y and the angle around the cylinder corresponds to the azimuthal angle the polar angular coordinate Q is 
suppressed. The blue dot is the base point x. The continuous red line is the null wavefront, emitted at time T = 0, at the time 
T = Tq. The dashed red line corresponds to the trajectories up to time T = Tq of representative null geodesies. The region 
shaded in green corresponds to the intersection of the normal neighbourhood of x with the hypersurface {T = To}. In (b), the 
region shaded in black at the bottom of the figure corresponds to spatial points which, at the time T — To, have been reached 
by the wavefront after crossing one caustic point. This region lies outside the maximal normal neighbourhood of the base point 



IV. EVALUATION OF Gr(x,x" 
A. The basic decomposition. 



We assume henceforth that ( 19 1 holds. Then, as noted, Gb.{x' , x") and Gii{x, x') can be written in Hadamard form: 

Gr{x\ x") = [C/i<5(ai) + V^H{-ai)]H{t' - t"), (26) 

and 

Gr{x, x') = [C/2<5(f72) + V2H{-a2)]H{t - t'), (27) 

where 

a, = a{x",x') = -lvl + ll! = -l{t'-t'r + l{y'-y"r + lll (28) 

^2 = a{x\x) = -\vl + lil = -lit-tT + l{y-yr + lii (29) 

with an obvious definition of 771,72, etc. Note that the subindices 1,2 no longer refer to the different 2-dimensional 
submanifolds as in Section III. We also have 



1/2 



A;=0 



Our aim is to calculate Gr{x,x") as given by ([t]). The following observations are of use in this calculation, and 
will be used frequently and without necessarily being flagged. 

(i) Gr{x,x") vanishes when x and x" are spacelike separated (i.e. when there is no causal curve from x" to x). 
Thus we are concerned only with points for which 

a{x",x)^^\{t-t"f + \{y-y"f<^- 
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Equality can be treated as a limiting case. In the (general) case of strict inequality, this indicates that the 
projections of x and x" are timelike separated in M2. Then by a Lorentz transformation in M2, we can arrange 
y" — y. This is assumed henceforth. 

(ii) Similarly, by a translation in M2, we can assume that y" ~ y ~ 0. 

(iii) By a rotation, we can take n{x") to be at the north pole, so that 9" = 0. {t^{x) is the projection of the spacetime 
point X onto §^.) Notice then that 71 = 0' . 

(iv) By uniqueness, the Kirchhoff representation ([t]) of Gr{x, x") is independent of the choice of t' e {t" , t). It will 
be convenient to impose an equal time split by setting t' ~ {t" +t)/2 so that 

t'-t" = t-t' = =: -. 

2 2 

(v) The time derivatives dt'Gii{x' , x") and dtiGR{x, x') include terms proportional to S{t' — t") and 6{t — t') respec- 
tively. But given the conditions ( [l9| , these distributions are identically zero for all spacetime points considered. 
Likewise, H{t' - t") = H{t -t') = l throughout. 



Applying these simplifications and using ([26]) and ([27| in we find 

6 

47r 



1 ^ 

G^[x,x") = ~—Y,Gk, (30) 



k=l 



where 



Gk ■■= / Gkdr, 1 < fc < 6 (31) 



and 

Gi ■■= lUiU2{S{ai)5'{a2) + S'{ai)S{a2)), (32) 

02 -^{UiV2 + U2Vi)5{cTi)Sia2), (33) 

^3 ^{UiV2d'{ai)H{-a2) + U2ViS'ia2)H{-ai)), (34) 

Gi := Ui{drV2)d{c7i)H{-a2)-U2idt'Vi)S{c72)H{-ai), (35) 

^5 -^ViV2{6iai)Hi-a2) + S{a2)H{-ai)), (36) 

Ge iVidt'V2-V2dt'Vi)H{~ai)Hi-a2). (37) 

These are grouped so that (in terms of the coefficients of the distributions) Gi is of the form "[/ x [/" , G2 — G4 
are of the form "[/ x V" and G5 and Gg are of the form x V" . In the following subsections, we evaluate the 
contribution of these six terms to Gr{x,x"). Figj2] illustrates the different contributions to Gr under the Kirchhoff 
integral method. 

Before proceeding with these calculations, we note the following regarding coordinates on the §^-component of 
S^/ = M X S^. It is convenient to use the geodesic distances of 7r{x') from Trix) and 7r(x") as coordinates on this 
So we introduce coordinates (a, /3) on given by 

a := 7i(^'>';e">") = e', 



Using (17l - with the appropriate replacements made for 72 - we can show that the invariant volume element on §^ is 

d2v'=Lo{a,/3;e)dadp, u{a, P;0) -.^ — — --^ — — , (38) 

(sm a — cos^ p + 2 cos cos a cos p — cos^ t^)^/^ 

and that for a function Q intcgrablc on S^, 
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FIG. 2: (Colour online). Diagram representing the M2 x §^ spacetime and illustrating the Kirchhoff integral method. The 
time coordinate T runs along the vertical axis, the angular coordinate 7 (and its negative, —7) runs along the horizontal axis 
and the Y coordinate runs perpendicular to the T-7 axes. Hypersurfaces at three different values of T are represented. The 
red curves correspond to the null wavefront ai = based at x" , while the orange curves correspond to the null wavefront from 
X. Their overlap on the hypersurface {T = t'} is shown with different colours corresponding to the different contributions to 
Gr under the Kirchhoff integral method ( |31[ ): red dots correspond to the 'f/ x U" term, the green curve (to the right of the 
dots) and the black curve (to the left of the dots) correspond to different "(7 x V" terms and the shaded region corresponds to 
the "V X V" terms. 



Qd2v' = 



da / Q{a, (3)uj{a, P)dl3 
/o Je-a 

da / Q{a,/3)uj{a,f3)d/3 
Ja-e 

da I Q{a, /3)uj{a, f3)df3 

iTT — ti Ja — 9 



(39) 



In fact this decomposition applies only when < < 7r/2: a similar decomposition applies when 7r/2 < 6 < n. 
The calculations below are done for the former case only, but have been checked to apply also to the latter case. 
We provide a separate (and much shorter) calculation to deal with the case 7(0;, x") = tt. This yields a significantly 
different result: see subsection F. 

For the upcoming calculations it is important to determine the singularities of the function uj{a,/3]9) in (38). This 
function is singular when 

cos^ = cos(a ± 6), 
that is, when /S = f3i + 2rm for l<i<4, nGZ where 



^4 



(40) 



We note that the function a;(a, /3; 9) diverges near these singularities like |/? — /3i|~ ' , i = 1, 2, 3, 4. 

Since the calculations are lengthy, we describe here the main features of the strategy that we follow. We seek to 
evaluate the triple integrals of plj ). The presence of delta functions will sometimes reduce the dimension to two or 
one. However, all the integrals we encounter share the feature that the integrand becomes infinite, typically at the 
boundary points of the domain of integration. The strategy we follow is to separate the integrand into singular and 
non-singular terms. For the latter, we will be able to show that these contribute overall finite terms to Gr. For the 
former, we introduce regularisation procedures that allow us to determine the distributional term corresponding to 
these singular integrals. 



10 



B. The U X U contribution. 



We have 



We recall that 



Gi= f GidV = f dy' f gid2v'. 



(41) 



f{x„)=0 



S{x - Xn) 



for any difFerentiable function /. Then, noting again that t' — t" > 0,t — t' > 0, we can write 

6{a,) = -{Siy' +y,)+d{y' -y,)}Hit' -t" -a), 
2/1 



6{a2) = -{Siy' + y2) + Siy' -y2)}Hit-t' -(5), 

2/2 



where 



2/1 

2/2 



(42) 
(43) 

(44) 
(45) 



Then we can carry out the y' integration in (41 ) - where we note that the integrands of the two terms concerned are 
both even functions of y' . We find 



dy' I U{a)UmSiai)S\a2) + S{<Ji)5'ia2)]d2v' 

S'i(T2)\y'=y, 



§2 



/S2 

U{a)U{P) 



2/1 



-H{t - 2a) 



2/2 



H{t - 2/3) 



d2v' 



= 4 



§2 



^(«)^(/3) 
a/? 



H{t - 2a) H{t - 2(3) 
- 4a2 - 4/32 



a45(/3-a)}rf2^^' 



§2 



5/3 I w(a,/3; 



F(a,/3; 0)9/3 



^ ap 

H{r-2P) \ 
Vr^ - 4/3' J 



H{T-2a) H{t-2(3) 
\/t2 - 4a2 ^r2 - 4/32 

i5(/3 — a)dadf3, 



8{I3 - a)dadp 



(46) 



J{a,P;e,T) 



where 



F{a,/3;d) := w(a,/3;l 



a/3 



Applying the decomposition (39), we can carry out the /3— integration to obtain 

'6/2 



Gi =4t 



9/2 



J{a, a; 6, T)da. 



(47) 



(48) 



We know that Gr{x",x) - not necessarily as generated by ([7|) - is a distribution, and so we should not expect the 
integral (48) to converge to a function. However we do expect it to be a distribution. This is indeed the case, and to 
identify this distribution, we introduce the following regularization procedure. From ( 40 ) we see that the integrand 
is singular at the limits of integration, therefore we let 



Gi = lim G ' 



(49) 
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where 



Integrating by parts yields 



where 



G^'^ := 4r / J{a,a;e,T)da. 

9/2+e 



•7r-e/2-e 



-Ji + J2-I 



(50) 



Ji 
J2 
I 



y/r^ - (6 + 2e)2 O^/sin^t 
H{t -2n + e + 2e) 



1 



^t2- (27r-6l-2e)2 (27r - 6')%/sh^ 
^-'/^-^ _Qia,9) 
9/2+t 



4^2 



H{t - 2a)da 



and 



Q{a,e):= —{Fia,a;e)}, F(a,a;( 



\/2sin(6'/2) a(cos e - cos 2a)i/2 ' 



We note that for each fixed 9, Q is analytic on (|,7r 



|), but is singular at the endpoints of this interval. 
(|5l|) and (|52|), O is the usual Lifschitz big-Oh, and the asymptotic relation holds in the limit e — > 0+. 



(51) 
(52) 
(53) 

(54) 
Also, in 



Removing the step function from the integrand of / yields 



/ = H{t ~e~ 2e)H{-T + 27r - 6* - 2e) 



Q{a,9) 



B/2+t 



4a2 



da + H{t -2TT + 9 + 2e) 



e/2^ 



Vr2 - 4a2 



da. (55) 



It is appropriate to pause at this juncture to consider the result that we expect to obtain for Gi. Given that Gi 
is constructed from the singular parts of Gii{x' ,x") (field) and Gr{x,x') (propagator), we expect it to correspond to 
the "most" singular part of Gr{x,x"). On general grounds, we expect that the spacetime location of the singularity 
in Gi is a subset of the future null cone of x" . That is, we expect that Gi{x,x") is singular at points x which are 
connected to x" by a null geodesic. Using the notation above, and recalling that 6" = 0, these are points for which 



= (-1) 



2n7r, n > 0. 



Here, n counts the number of caustic points through which the geodesic has passed: for n = 1, the geodesic has 
passed through the caustic x with Tt{x) antipodal to 7r(a;"); for n = 2, it has passed through this caustic and that 



with tt{x) = tt{x"), etc. From the form of G'f'^ above, and the expression (55 1 for /, it is clear that these are 'points 



of interest' in our investigation. We consider the case of most interest, which corresponds to r e (tt, 27r). In this case, 
X lies outside the maximal normal neighbourhood of x", and the null geodesic spray from x" has formed exactly one 
caustic. Then n — \, and we expect to see that Gr{x,x") is singular when, and only when, t = 2Tr — 9. 

We recall that our aim is to calculate the singular part of Gu{x, x"). In this calculation, we will frequently encounter 
terms A that are finite for all values of the spacetime variables (perhaps of the region of spacetime presently under 
consideration, e.g. tt < t < 27r), and for all sufficiently small values of e > 0, and that remain finite in the limit 
e — >■ 0+. Such terms do not, by definition, contribute to the singular part of Gr- We use the following notation to 
characterise quantities which possess these features: 



A = 0{1). 



Using (55) and noting that t > t: > 9, and introducing (for ease of notation) a :~ 2a, we have 

/ = H{-T + 2n-9- 2e)Ii + H(t - 2tt + 9 + 2e)l2 

where 



(56) 
(57) 

(58) 
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and 



R{a, 



V2sm{e/2) 



(acos(a/2) — 2 sin(a/2))(cos — cosa) — asin(a/2) sina 
a2(cos0 — cosa)3/2 



(59) 



Notice then that R{a, 6) is singular at the points a — 9 and a — 2tt ~ 6. These correspond to the lower and upper 
limits of I2, and to the lower limit of /i, when e ^ 0+. Further, a = 2Tr ~ 6 also corresponds to the upper limit of Ii 
in the limiting case r = 27r — 0, e — >■ 0"*". So in order to capture the dominant contribution to /, we need to isolate the 
singular behaviour in R{a, 6), leaving a finite remainder that will contribute an overall 0(1) term. In fact we proceed 
by isolating the singular behaviour of R{a, 6) /\/t + a. This greatly simplifies the resulting integrals. Notice that the 
factor is analytic on [9, 27r — 6] and so does not introduce any additional singular behaviour. 

We define 



lim (a 

-^6 + 



lim (a 

a^0 + 



Rli9,T) 
R2{0,t) 

R3{0,t) 
Ri{9,T) 

Of these, the following values are required 

Ri{e,T) 

R3{9,t) 
Ri{0,T) 



^3/2 



R{a,9) 
y/r + a' 
R{a,9) Ri{9,t) 



R{a,9) 



lim {2n-e-af/^ ^ , 



lim {2Tr-e- a)i/2 

a^(27r-e)- 



1 

V29Vsh[9V7T9' 
1 



i?(a,( 



R3{9,t) 



y/r + a {2tt 



a)3/2 



y2(27r - 9)Vsh[9y/T + 2tt - 9 
{T + 2n- 9) (27r - 6*) (2 + cos 9) + 2(2r + 27r - 9) sin t 



4^/2(27r - 9Y{t + 2ti- 0)3/2(sin6')3/2 



(60) 
(61) 
(62) 
(63) 

(64) 
(65) 
(66) 



The Ri are the coefficients of the divergent parts of Rj \Jt + a at the different endpoints. It follows that we can write 
R{a,9) Ri{9,t) , R2{9,t) , Rz{9,t) , Ri{9,T) 



y/7T^ {a-9f/^ (a-6')i/2 {2Tr - 9 - a)^/^ {2n - 9 - a)^/ 



— +i?(a,0,r). 



(67) 



where for each r S (tt, 2tt) and each 9 G (0, tt/2), R{a, 9, t) is a continuous function of a on the interval [9, 2tt — 9\. It 
follows that 



^ R{a,9,T) 
\Jt - a 



da — 



6l+2e 



2ir-e-2e 



e+2e 



R{a,9,T] 
\Jt - a 



da = 0(1). 



(68) 



That is, R makes only an 0(1) contribution to /. The following antiderivatives are then needed (constants of 
integration are omitted): 



da 



2^/7^ 



V^^ia - 6')3/2 (61 - T)Va^' 
da 2\Jt — a 



da 



V7^{a-9y/^ 



= 2 arctan 



da 



\Jt - a 



(69) 



V7^(27r - 61 - a)i/2 



V7^(27r - 61 - a)3/2 {r - 2n ^ 9)^2-!: - 9 - a 
Applying these (and using ([68|), we see that all terms proportional to Ri are 0(1), and we obtain 



= -2 In WV^^ ^j2-n-9-{Sf) 



- T)^/2t (r - 27r + 9y2-K - 29 - 2e 

2R4, |ln(\/27r -9-t)- ln(Vr - 9 - 2e + V2n - 26* - 2e)| + 0(1) 



(71) 
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and 



h = -2i?i 



^/T-6-2t 



2R. 



t)V2^ (T-27r + i 



Vt - 277 + 6» + 2e VT-0-2e 



/2e 



V27r - 26* - 2e 



-2i? 



4 |ln(%/T - 27r + 61 + 2e + \/2i) - ln(Vr - 61 - 2e + V27r - 26* - 2e) | + 0(1) 



(72) 



Now we collect terms from (151), (152|), (|7ll) and M evaluate g['^ via (jSOl) and (|57|). 
We note first that the Ri term in / cancels to 0{e) with Ji (using t > 9 in Ji), and so there is only an 0(1) 
contribution to Gi from these terms. Noting that the i?i contribution to / arises by formally setting R2 — Rj, — Ri — 
0, we can write 



lim (Ji + /|fl,=i?3=i?,=o) - 0{l). 

e— >0+ 



(73) 



The contributions from J2 and the i?3 terms in / arise by formally setting i?i = i?2 = i?4 = in /. These yield 



1 

(27r - e)yJsme^/T + 27r 
+ iJ(r - 27r + 6* + 2e 
+0(1) 



X <^ i/(-T + 27r - 6* - 2e) 



Vt - 6* - 2e 
Vtt - - e 



1 

T-27r- 



Vr - 27r + I 



2e 



Vr - 27r + 6* + 2e 
T - 27r + 6* 



1 Vt - 61 - 2e 

h 

v/e Vtt - 6* - e 



r-27r- 



To evaluate the terms in braces here it is convenient to introduce 2; := r — 27r + 0. Note then 

Vt - 61 - 2e 1 \/27r - 261 - 2e + z 1 



r - 27r + I 



Vtt - 61 - 
+ 0(1), 



(74) 



and so we can write 



where 



(27r - 0)Vsh[eVT + 2TT-0 



A(^)(z)+0(1) 



aW(z) :=^i/(-z-2e) + ^(l-^^ 
^ z ^ z \ Vz + 2e 



H{z + 2e). 



In Appendix A.l, we derive the distributional identity 



lim aJ''(z) = V2PV ( - 



e->-0+ 



where PV is the principal value distribution, defined by the following action on a test function 



< PV 



>= lim 

£->0 + 



(^) 



dz. 



Thus 



v/2 



1 



^nm (J2 - /|«,.«.=«,=o. _ e)V^eV r + 2n-e ^ ' \r-2n 



0(1) 



(75) 



(76) 



(77) 



(78) 



Finally, we consider the i?4 contributions to G\ . These take the form of logarithms. We note that in (71) and 



( 72 1 , due to the range of values of r and 9 under consideration, we have 

In \Vt-9-2€ + V2tt -29- 2e\ = 0(1). 
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We also have the distributional identity (see Appendix A. 2) 



It follows that 



lim ln(Vz + 2e + V2e)H{z + 2e) = \n{^/i)H{z) 

e-i-0+ 



lim I\r,^r^^r,=o = -i?4 In |r - 27r + 6l|(i7(-T + 27r- 61) + i?(r- 277 + 6*)) 

e— i.0+ 



-i?4 ln|r-27r- 



Collecting the results of (73), ( |78[ ) and (801 completes the calculation of Gi: 

Gi = hm G\ = , X 

£^0+ (27r - 6')%/sh^V'r + 27r - 6* 



PV 



1 



r - 27r + 6* 



(r + 27r - 9) (27r - 61) (2 + cos 61) + 2(2r + 27r - 9) sin 61 



(79) 



(80) 



%/sin6lV27r- 61 



PV 



T-2n + e 



8{2tt-9){t + 2tt-9) smO 
(27r- 61) (2 + cos 61) + 3 sin 61 
8(27r-6l) sin6l 



In It -2n 



ln|T-27r + 6l| } + 0{1) 



In the last line here, we have used f{x)PV{^) ^ f{0)PV{^) + 0{1) and similar with PV replaced by In, 
differentiable function /. That is, we set r = 27r — in the coefficients of PV and In. 



(81) 
for a 



C. The V xV contribution. 



Two terms contribute to Gr(x,x") under this heading, corresponding to G^ and Gq. We show here that 

G5 = 0(l), G6 = 0(l) (82) 



for all T G (0,27r). 



We consider Gg first. Using the functional forms of Ui and Vi [i = 1,2) given in (281, (29) and (22), we see that 

T 



dt'Vi - - — 9,/ 14, dt'V2 = —dy'V2. 
2y' 2y' 



(83) 



Then combining (31) and (37) gives 



Gfi 



Noting that each Vi is an even smooth function of y', we see that this is the integral of a smooth function over the 
compact set {x' e Sf : cri{x,x') < and a2{x',x") < 0}. Hence 



G6 = 0(l). 



(84) 



Prom (31) and (36), and using (42) and (431, we have 



Gs - -T / d2v' r dy' 

JS2 Jo 



^<5(2/' - y,)H{T - 2a)H{~a2) + —5{y' ~ y2)H{T - 2P)H{~<t,) 

yi 2/2 



-2r / d2v' 

JS2 



Vt^ - 4a2 



v/t2 - 4/32 



where yi and y2 are given in (44| and (45). The terms i^o{a)V2\y'=y-i and iyoif^)^i\y'=y2 analytic on S , and so both 
integrals are finite: the integrand diverges at r = 2a and at r = 2/3, but only at the rate (r — 2a)^^^^ (respectively 
(t — 2/3)~^/'^), yielding an overall integrable term. Hence we can state that 

G5 = 0{1). (85) 
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D. The U X V contribution. 



The terms G2, G3 and G4 form the U x V contribution. From (28) and (29), we see that 



Using this and ( 83 ) , we can combine ( 34 ) and (351 to obtain 



G3 + G4 = 



§2 



d2v' / dy' 



^^^H{-a2)dy, {V25{a,)} + 1 o 2 



(86) 



where here and below, 1 -s-)- 2 indicates dupHcation of the preceding expression with the replacements 1 o 2 and a o /?. 
The pole at the origin shows that, as above, the integral must be regularized. We write G2 + G3 + G4 = \\m^^Q+ g'^\ 4, 
where 



G 



3,4 



d2v' 



§2 



dy' 



^H{-a2)dy, {V2Sia,)} + 1^2 



(87) 



To proceed, we integrate by parts. This introduces three types of term: boundary terms corresponding to y' — ±e 
(these are collected in Ki below) , terms arising from the derivative dy' acting on 1 /y' (these are collected in K2 below) 
and terms arising from the derivative acting on H{—fJi). These last yield terms with a product of delta functions in 

the integrand. These terms cancel identically with the Gj^'' contribution (with the obvious definition of g'"2^ - the 
integral (31 ) with k = 2 and a principal value truncation on the y' axis). Thus we can write 



G2 



where 



with 



Ki := -T / d2v' 



§2 



dW 



U{a) 
U{a) 



G3 + Gi= lim Gi% 4 

e— i.0+ ' ' 



Gi%^ = K,+K2, 



H 



= -2t 



d2v'U{a) 



V2 






' t2 


£2 


a" 


e 


y 


V 


'2- 




V2 




H 


' t2 






e 


y' 


V 


'2- 




V2 




A 




£2 




e 


y' 




-2' 


2 



2 



1 o 2 



and 



K2 



d2v' 



§2 



d2v' 



§2 



= 2t 



doV 



H{a- (i) + 1 o 2 



dy' y^^^5{cr^)H{-a2) + 1^2 

L y' 

^^V2\y'=y,H{T - 2a,)H{a - /3) + 1 o 2 
,U{a) 



(89) 



V2\y'^y,HiT-2a,)H{a- (3), 



with 



(90) 



(91) 



The first expressions for each of Ki, K2 follow by noting that the integrand is an even function of y' . The second 
expression for Ki relies on the fact that there is a contribution only at the zero of the argument of the delta function, 
and that a,li are both positive so that H{a^ — /S^) = H{a — /?). The final expressions arise because the 1 ■(-> 2 
contributions are equal. (This last fact is not immediately obvious, but follows from a decomposition equivalent to 
(39) in which the order of integration with respect to a,/3 is reversed.) 
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It is useful to consider the general form of the integrals Ki,K2: using the decomposition (39), it can be shown that 
for any function / that is integrable on 



d2v'f{a,f3)H{a- /3) 



+2 
+2 



da 

e/2 Je-a 

■rr-e/2 

da 



fia,P)ujia,p)dp 

fia,l3)Lu{a,P)dl3 

9 



da 

TT-e/2 Ja-6 



f{a,(3)oj{a,f3)dp. 



We will then naturally encounter these terms: 

Li(a) 
L2{a) := 
Lsia) := 
where 

oo 
n=0 



uj{a, /3)i'{a, I3)dl3, 

— a 
a 

uj{a, f3)i'{a, I3)dl3, 



uj{a, l3)h'{a, P)dj3, 



The term Ki corresponds to an integral of the form ( 92 ) with 



v. 



4t 



U{Tj2)6{a-Tj2)V2\ 



y' —e,a.—Tf /2 



4t. 



= U{Tj2)S{a-Tj2)v{Tj2,P), 



where 



Then using (92) and the definitions (|93|)-(|95|), we find 



8r 



C/(Te/2) {i7(r, - - r,/2)Li(r,/2) 



+iJ(r,/2 - 0)i/(27r - - t,)L2{t,/2) 
+ iJ(r, - 27r + 6l)i7(27r - t,)L^{tJ2)} . 



Likewise, 



K2 = 32t 



(7(a) 



,/2 (r2 - 4a2)3/2 

e (t2 - 4a2)3/2 
C/(a) 



Li{a)H(T - 2a^)da 
L2{a)H{T ~ 2a,)da 



.-9/2 {r' - 4a2)3/2 



L3{a)H{T - 2a,)da 



(92) 

(93) 
(94) 
(95) 

(96) 



(97) 



(98) 



(99) 



(100) 
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The terms H{t — 2af_) here arise naturally, but for the next steps we note that H{t — 2a^) = H{t^/2 — a). This 
follows from the definitions (91 1, (98). The step functions can be extracted from the integrals to obtain 



K2 = 32r <^ H{t, - e)H{e - re/2) 



U{a) 



J2 



U(a) 



(r2 - 4a2)3/ 



-Li{a)da 



U{a) 



-L2{a)da 



H{n - 27r) / 

J 7T — 



(r2 - 4^2)3/2 

■^-'/^ U{a 

(r2 - 4a2)3/2 

U{a) 
e/2 (t2 - 4a2)3/ 



L2{a)dc 



+^(--/^-^)i M^--4;2)3/2 ^^M^" 

+H{t, -2Tr + 0) f 
Je 

rrc/2 

+H{t, -271 + 0)H{2tt - T,) / , , ^Y'., L3ia)da 

JTT-e/2 [t^ - ''-■2^0/2 



_e,2 {r^ - 4^2)3/: 



■;L^{a)da 



(101) 



We consider from now on the the case of most interest, when tt < t < 27r, so that the point x lies outside the 
maximal normal neighbourhood of x", after the formation of the first caustic but prior to the formation of the second 
caustic. Hence for sufficiently small e, we have n < < 2n and since < f? < 7r/2, we also have > 26 > 9. Then 
some terms of (99) and (101) are identically zero and others simphfy, leaving 



where 



and 



with 



8t 
8t 



U{Tj2)Hi2Tr - e - n)L2{Tj2), 
U{Tj2)Hin-27T + 0)L3{Tj2) 



K. 



2{i) 



K. 



2{ii) 



K2(iii) + K2(iy) 



(102) 

(103) 
(104) 

(105) 



K. 



2{i) 



K. 



2{ii) 



K. 



2(iii) 



K. 



2(iv) 



^^^L/2 (r2 ^2)3/2 



tJ2 



32tH{2tt-9-t,) 

Je 

S2tH{t, -2tt + 



U{a) 



g (r2 - 4^2)3/2 

iT-e/2 



U{a) 

'g (r2 - 4^2)3/2 

U{a) 

^^g/2 (t2 - 4a2)3/2 



L2{a)da, 
L2{a)da, 
L^{a)da. 



(106) 
(107) 
(108) 
(109) 



The nature of the integrals contributing to Ki and K2 depends on whether or not the singular points ( 40 ) of 



the function w(a, (3) occur in the domain of the integrals. In some of the cases below, we see that the singularities 
occur only in certain limits. It is nevertheless crucial to take account of these, as the overall singular contributions 
to Gr{x,x") will include contributions from these limiting singular values. As with the calculation of the U x U 
contribution, we follow the strategy of decomposing integrands into singular and non-singular terms, which lead 
respectively to singular and non-singular 0(1) contributions to Gr. 
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1. Evaluation of K2(i)- 



This term involves Li{a), given in (93|. We note that z/(a,/3), defined in (96l, is an analytic function of (a,/3) for 



all (a,/3) in the relevant range. The integrand of Li is singular only at the lower limit, with 

uj{a,(3)^O{{P-0 + a)-^/^), 13 ^ {d - a)+ . 

Thus Li{a) is finite for all a £ [^,0]. Furthermore, t > 2a for all a G [^,9], and so -K'2(i) corresponds to the integral 
of a continuous function over a finite interval, yielding 

X2W=0(1)- (110) 

2. Evaluation of Ki(^ii-j + K2(ii) ■ 



We have 



«/2 



U{a) 



g (r2 - 4a2)3/ 

with L2 as given in ( [94| ), and note that a e [0,r(:/2]. Then oj is singular at the lower limit (3 — (5^ of the integral L2, 
and is also singular at the upper limit /? = a in the limiting case when a = rg/2 and = 2tt~6. No other singularities 
arise in the integral. To take account of both singularities, which correspond io (3 = a — 6 and (3 = 271 — a — 9, we 
define 



sinasin(a 



1/2 



gi(a):= lim (/3 - a + g)V^c.(a,/3) = ( — — ) (m) 



.2(a) := lim , (2. - a - ^ - ^) ^^a, ^) ^ f -^"t"^"/ V'' . (112) 



and 



,g-i-(27r-a-e)- \ 2sin0 

Then we can write 

^a,/3) = + ^=J^M_ +a.2(a,/3), (113) 

VP — a + p V27r — a — — p 

where, for all a in the relevant range, ^2(0?, /3) is a continuous function of /?, and daUJ2{oi^ (3) is an integrable function 
of j3. By expanding the analytic function v{a,/3) about the relevant end-point, we can then write 

L2ia) = r \^M=[uia,a-e) + Oi(3-a + 9)] 

^^^"^ 'iy{a,2TT-a-9) + Oi2Tr-a-9- l3)]+u}2ia,P)iy{a,P)\dp 



y/2TT -a- 9- 13 

dp 



qi{a)i'{a, a — 6) 



a-e V/3 - a + I 



+q2{a)v{a,2TT - a - 6) \ ^^^+£2.2 (a) 

y/2'K - a - - j3 



= 2V9qi{a)i^{a, a - 6) + 2^^^q2{a)iy{a, 2tt - a - 9){Vt^^ - v^tt - 9/2 - a)) + i2,2(a), (114) 
with an obvious definition of a continuously differentiable function £2 2- It follows that 



L2{a) L2{a) + 2^/^q2{a)i^{a, 2tt - a - 9)y/TT~9/2-a (115) 
is on [9, T^/2] for all < 2tt — 9. We can then calculate 
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2{ii) 



"-/^ 1 U{a) 

e (r - 2a)3/2 (r + 2a)3/2 

1 



r,/2 



, (r - 2a)3/2 

_23/2 / 



e (r + 2a)3/2 



(2t)3/2 



L2{a) - 2^l^q2{a)v[a, 27r - a - 0)^/1: -9/2^ a 
da 



^^^/^^Z2(T/2)+0(r-2a) 



(^^2(r/2) 
_23/2t^(^/2) 



da 



(2t)3/2 

U{t/2) 



(2t)3/2 



, (r - 2a)3/2 
q2{T/2)y{T/2,2^-T/2- 
1 1 



y/TT - 0/2 - . 

^ (T-2a)3/2 

'/^ V^-0/2-< 
(r - 2a)3/2 



da + 0(1) 



-da + 0{l) 



L2{t/2) 



^/^^^c - 29 

2^/\2(T/2)j/(r/2, 27r - t/2 - 



y/TT-9/2-Tj2 1 



ln|r-27r + 6l| 



23/2 

^L2{T/2)^J7/2e-^ + q2(T/2)iy(T/2, 27r - t/2 - 61) In |t - 27r + 6l|| + 0(1) 

Incorporating the expression for A'x(ii) from (103), we find that 

.U{r/2) 



0{l) 



(2t)3/2 



lim (ifi(ji) + if2(ji)) = 16- 



See ( 22 1 , ( 23 1 and ( 96 ) for this last equation. 



3. Evaluation of K2 (Hi) ■ 



The calculation of this term is very similar to that of K2[u). We have 



K. 



2{iii) 



32tH{t^ - 2n + 9)K2(iu), K2(iii) 



r-e/2 



U{a) 



g (r2 - 4a2)3/2 



L2{a)da 



(116) 



-q2{T/2)v{T/2, 27r - t/2 - 9) In |t - 27r + 9\H{2t: - 9 - t) + 0(1) 

= 16 ^ /'""^^^l g2(7r - 6i/2)i/o(7r - 61/2) In |t - 27r + 9\H{2tt - 9 - t) + 0(1). (117) 
V2(27r - 9) 

We note the use here of /(t) In \t - 2n + 9\ = /(27r - 0) In \t - 2n + 9\ + 0(1) for a function /, and of 

via, a) — voia). 



with L2 given by (94). We encounter singularities of uj at the lower endpoint (3 — a — 9 and at the upper end point 

on c 

d/S 



/3 = a in the limiting case a = t: — 9/2. Thus we carry out the double end-point expansion of (113) and write 

d/3 



-^2(0:) = qi{a)v{a,a 



Q-e - a + ' 



+ q2{a)v{a, 2n — a 



_g V27r - a 



+ L2.2{ol) 



= 2V9qi{a)iy{a, a - 9) + 2^/'^q2{a)v{a, 2ti - a - 9){Vn - a - \J n - 9/2 - a) + ^2,2(0), 

where L2,2 is on the domain of L2. It follows that 

L2{a) = L2(a) + 2^/^q2{a)v{a, 2tt ~ a ~ 9)^/tt -9/2 -a 

is also on this domain. Thus we can evaluate K2{iii) by expanding UL2/{t + 2q;)3/2 around the upper end-point 
of the integral - we note that this guarantees that we capture the most singular and hence dominant contribution to 
the integral. (Note also the slight difference at this point in comparison to the calculation of K2{ii)-) 

U{a)L2{a) _ U{t:-9/2) 



(T-K2a)3/2 (r + 27r - 61)3/2 



|L2(7r - 61/2) - 2^/2g2(7r - 9/2)iy{7r - 61/2, tt - 61/2) ^tt - 61/2 - a + 0(7r - 61/2 - a)| . 
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Hence 



~ t/(7r-g/2) f"'^' 

0(^-0/2 -a) , 



(r + 27r - e')3/2 ' 'WT-27r + e' ^f^^ 

U{t: - 0/2) 



{t + 2tt~ 6')3/2 



q2in - e/2)va{Ti - 6/2) In |t - 27r + i 



+0((^ _ 2^ + 0)1/2) + o (^^^^) ■ (118) 

Note the use of v{a, a) — vo{oi)- Given that r > tt > 20, we can therefore write 

lQU{TT-e/2)L2(T^~e/2) 

+ """'",'^"'7''"' -°''' - 0/2) I- - 2' + "i^f- - 2' + ") + -^d)^ (119) 



4- Evaluation of Kn^i^-^ + A"; 



2(iti) • 



We have 



i^2(^.) = 32TiJ(T, - 27r + 0)if2 



K. 



2{iv) 



./2 



t/(a) 



._e/2 (t2 - 4^2)3/: 



-L3(a)iia, 



with I/3 as defined in ( 95 ) . The integrand of is singular at the endpoints of the interval of integration, which, in 
terms of the definitions above, are (3^ = a ~ 6 and f3i + 2tt = 2-k — a — . We note that ^2 < for the range of 
values of a of concern here. Also, with a > tt — 0/2, we have a + > tt + 0/2 > tt, and so27r — a — 0<a + = P^. 
Therefore /Ja and /3i + 27r are the only points at which the integrand of L3 is singular. So again we use the double 



end-point expansion (113) and write 



L,{a) = r " ' \^fi^^Ha,a-e) + 0il3-a + e)] 
\-e L V P - a + fc* 

-=J^M== [v{a, 27r - a - 0) + 0{2^ - a - - /?)] + W2(a, /3)z^(a, /3) |> d/3 
V27r — a — 6/ — p 



(7i(a)j^(a, a — 0) 



/'27r — Q — 

-|-g2(a)i^(Q!, 27r — a 



dp 



9 V27r - a - - /3 
-1-^3,2 (a) 

= 23/2qi(a)j/(a, a - 0)(7r - a)^/^ 2^l'^q2{a)v{a, 27r - a - 0)(7r - a)^/^ 2(a), (120) 

with an obvious definition of L3.2 By analyticity of v and the properties of uj2 noted above, it follows that L3.2 is a 
continuously differentiable function of a on [tt — 0/2, r/2]. 
Then we calculate -ft'2(ii,) by expanding -^3(0;) about a — r/2: 

L3(a)=i3(T/2) + 0(a-r/2), 



21 



so that 



2(iv) 



f Uir/2)L,ir/2) 
'^_e/2l(2r)3/2(r-2a)3/2 

U{t/2)L3{t/2) r 1 _ 1 
(2t)3/2 [V^^e V^^2^ 



+ 0((T-2a)-i/2)|da 
hO(l). 



(121) 



Feeding through to K2(iy), we see that the first term here cancels (to 0(1)) with Kn^^y-^ in the limit e — > 0+, and 
we find 



lim (Js:i(„) + K2(ty)) 



IW{t/2)L^{t/2) 
\/2rVT -211 + 9 
16C/(7r-6l/2)L3(7r-6l/2) 



H{t -2n + e) + 0{l) 

H{t -27^ + 6) + 0{1) 



^2(27r-6l)VT-27r + 

We use here /(t)/Vt - 27r + 6* = /(27r - 6l)/Vr - 27r + 6* + 0(1) for a function at 27r - 6* 



(122) 



5. Summary: the U x V contribution. 



Combining the results ( |110| , ( |117| ), ( |119| ) and ( |122[ ) and using 

L2{n-9/2) = L3{Tr-0/2) 



which can be seen to follow immediately from the definitions ( 94 1 and ( 95 ) , we obtain 



+ G3 + = m^L^lIM^Mn - 0/2) m \r-2. + 0\+ 0(1). 
^2(27r - 0) 

E. Gr{x,x") beyond the normal neighbourhood. 



(123) 



Combining (30), (811, (84), (85) and (1231 yields the expression for Gji(x,x") in the case where tt < r < 27r. This 
corresponds to the situation where x is outside the maximal normal neighbourhood of x" , and where the S'^-envelope 
of null geodesies emerging from x" have passed through a single caustic. We have 



Gr{x,x") 



Vsin 6*^/2 TT 
In \t-2tt + 1 



PV 
16 



r - 27r 



[/(vr - e/2)q2{^ - e/2)uo{^ - e/2) - 



(27r-e')(2 + cos6l)+3sin( 



47rV27r - 9 \ ^2 

Using the definitions (21), (23) and ( |112 ) of [/, vq and q2 respectively, this yields 

1 / 1 



2(27r6l)(sin(6l)3/2 



o(tm) 



GK(x,a;") = - 



7r\/sin6l\/27r - B 



PV 



r- 27r 



/27r-6 




1 — rn? 






V sine 







cot(27r - 6*) 



(27r-6l)2 



2n-e 



In |r - 2Tr - 



+ o(i) 



271-6* 
sin(27r - 9) 



1/2 



1 



We have used here linearity of PV, the property that /(0)PV^(l/a;) = /(a;)PF(l/.T) + 0(1) and the definition ([23|) 
of fQ. We have also absorbed an 0(1) term into the logarithm to facilitate the final step. This step is to rewrite the 
last equation in covariant form. We find 



Gnix, x") = U(2^ ^)PV ( / „. ) ~ i.o(27r - ^ ^^^^y 

\TTa{x, X") J TT 



(126) 
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F. The case 7 = tt. 



A perplexing feature of the result ( 126 1 is that it diverges at 7 = tt at points x not connected to the base point 
x" by a null geodesic. This can be understood by considering the method by which this result was obtained, and, 
in particular, the use of the property / = 0{1). Among the various terms swept under the 0(1) carpet are terms 
involving negative powers of sin0 (i.e. negative powers of sin 7). It is only permissible to ignore these when sin 7 is 
bounded away from zero. Thus the result (126) should be understood with this caveat in mind. 

Furthermore, the use of the coordinates (a, /3) on §^ is not valid when 7 = tt. This case corresponds to the use 
of the geodesic distances of a point on from the north and south poles as coordinates on S^, which is clearly 
invalid. However we can calculate G'fl(x, x")|.^=^ by returning to the calculation above prior to the introduction of 
these coordinates. To that end, we note that 



and so 



71 = a, 



1 2 , 1 /2 1 2 
8 2-^ 2 ' 



72 = /3 = TT - a, 



1 9 1/2 1 / x9 



As coordinates on we use the usual coordinates 9' = a and (j)' . The calculation of Gr{x,x") is made considerably 
easier by virtue of the fact that all integrands encountered are independent of (j)' . The calculation below assumes (as 
above) that tt < r < 27r. 

To determine Gi in the case 7 = tt, we take up the calculation leading to (46) at the second line: 

1/2 



Gl = T 



1/2 



§2 



= 2ttt 
We have 



(— ) 
V sm a / 

\/ a{TT — a) 



TT — a 
sin(7r — a) 



V' 



H{t - 2a) 



y'=yi 



y' 



H{t -2{tt - a)) 



y'=V2 



sin a da i 



5' (^2) 



H{t~ 2a) 



S'ia,) 



S'icji)\y>=y^ = ^^{(5(a-7r/2)}, 



iJ(T - 2(7r - a)) 



1 d 



da 



(127) 



S'i'y2)\y'=y,^~^^{S{a~7T/2)}. 



Using these and integrating by parts gives 
At d \ yJa{-K - a) 



Gl = 



TT Jq da 1 
At d \ yj a(7r — a) 

8t 



4a2 



H{t - 2a) 



y/a{Tr-a) 



\/t'^ — 4(7r — a)' 



.H{t - 2{tt - a)) ) S{a - 7r/2)da 



4a2 



H{t - 2a) 



ya(7r- a) 
y^r^ — 4(7r — a)' 



.H{t -2{tt - a)) 



Sttt 



= - (^2 _ ^2)1/2 '^(^ - ^) + (^2 _ ^2)3/2 - 

Next, we recall that 

G2 + G3 + G4 = lim G^ 

e-i.O+ 



(128) 



as given by (88). From (89), with 7 = tt (giving /? = a — tt and iJ(a — /3) = 1), we have 



-2t 
87rr 



,1/2 

e 



(— ) 
V sm a / 

r(asina)i/2 
7o e 



t2 



a 



sin a da ( 



Msin(r,/2) 



1/2 



y'=e 
00 

E 

fc=0 



(5(a - T,/2)da 
t.fe(7r-r,/2)(|(7r-T,))' 



(129) 
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Likewise, from (90), we obtain 
K2 = 



327rr 



(asm a 



[ (^2 _ 4^2)3/2 y2\y'=y,HiTj2 - a) da 
32^^ ( J!! 4a2)3/2 ^^''"^^ - «)(|(^ " da 



(130) 



We expand 



Q{a) — (a sin a) "'"^^ i^fc(7r — a) (7^(71' ^ 2q!) 



fc=0 

around a — t/2 and integrate. The leading order term is proportional to e^^, and cancels identically with Ki. The 
remaining terms are finite, and so 

Ki+K2 = 0{l), 

leading to a finite value for G2 + G3 + G4 in the case when 7 = tt. 

The analysis of Gs + Gg done in subsection C carries over to the case 7 = tt, and so 



G5 + G6 = 0(1) 



also holds in this case. 

Combining these results, we find 



2t 
TT 1 dr 



H{t - tt) 



(r2 -7r2)i/2 



H{t-tt) 



(r-7r)i/2(r + 7r)3/2 



+ 0(1), (r<27r). 



(131) 



Note that the final form here shows explicitly the distributional nature of the expression, which is not manifest in 
the previous line. Of course for r < tt, the 0(1) term is identically zero. 

This result shows that the only singularity arising at 7 = tt is confined to the null cone (as expected) . This highlights 
a drawback of the result (126): it does not apply on a neighbourhood of 7(2;, x") = tt. 



V. SINGULARITY STRUCTURE VIA MODE SUM APPROXIMATION 



In this section, we present the results of a completely different approach to the calculation of the singular part of the 
retarded Green's function outside the normal neighbourhood. This approach - involving asymptotic approximations 
for the special functions that arise in a mode sum decomposition of Gr - has the advantage of yielding global results. 
That is, we obtain the singular part of Gr globally. The drawbacks are that we must restrict to a particular value of 
the coupling term + ^i?, and we obtain results which have unresolved causality issues. (It is clear however that 
these issues must be resolvable, and they do not affect the correct identification of the singular part of Gr.) The 
results are of interest in that they are global in nature: they are not restricted either to the normal neighbourhood 
or to the region before the formation of the second caustic. What is of most interest is that they return the four-fold 
singularity structure discussed above. 

To begin, we find the Green's function solutions of the wave equation ^ by carrying out a multipole decomposition 
of a solution of the wave equation (|4]) . Imposing the appropriate boundary conditions we find that the Feynman [5J [TS] 
and retarded Green's functions can respectively be expressed as 

00 

Gf{x,x') = lim - V(^ + l/2)i/^''(A?7,)P,(cos7), (132) 

and 

00 

Gr{x,x') = 2Rc{Gf{x,x'))H+{x,x') = H+{x,x')'^{e+l/2)Jo (At;) (cos 7). (133) 
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The expression (132 1 is obtained by expanding Gp in terms of Legendre polynomials, and applying equation 
(2.77) of fT5l to obtain the resulting (l+l)-dimensional Feynman Green's function. In these expressions, rj^ := 

y i\t - f \ - itf -{y- y'Y and l{l + 1) + 2^. 

We find it more useful to carry out the large-^ asymptotic analysis on the Feynman Green function and then obtain 
from it the large-^ analysis for the retarded Green function. Therefore, we expand Eq.(132| for large-£ - i.e. large 
argument for the Hankel function ' - order-by-order. For that we make use of Eq.8. 451(4) of [17): 



1/2 



„-j(z-7r/4) 



n-l 

E 

.fe=0 



Oik 



1 r(fc-n/2) 

{2iY fc!r(-A: + 1/2) 



where ([l7]) P- 963) |02| < 1 if — 37r/2 < arg(z) < 7r/2 and Im(z) < 0, which is the case here: 



1\ At 

— e 
V 



e -> 0+. 



(134) 



lm(z) < then arises as we are only concerned with the case t] > (and note that At = \t — t'\). We have a 
contribution to Gr only when cr < 0, so that rj > 0. But ry = with cr < can only arise when x = x' , which case 
can be ignored. 

We take n — 2, giving 



1/2 



,-i(z-7r/4) 



1 



27 



Hz 



(135) 



Consider the remainder term GJ?™ - i.e. the contribution of 62 to Gp- This yields 



e^0+ 2^ 



1=0 



3/2 



'''/'^^2K+^) 6-'^"= (cos 7) 



(136) 



Since 16*21 < 1, we have 



00 

E^2(^+;.)-'/'e'^^''«P,(cos7) 



1=0 



<£(£-f 1/2)- 



■3/2 



1=0 



so that the series is absolutely convergent. Furthermore, each term in the series is continuous on = {(2,7) : \z\ < 
1,7 G [0,7r]}, and so the Weierstrass M— test yields uniform convergence to a continuous function. Therefore G^p™ is 
continuous for 77^ ^ 0. It follows that the singularities of Gf and Gn are contained in the terms 



G^p'^ := - lim 

^ e^0+ 2 



G 



(1) 



lim -W-e-/V^/'e-''«/25i, 

e-S.O+ 16 V TT 



(137) 
(138) 



where 



Sk :=^(^+l/2)i/2-'=e-'^''^Pf(cos7), = 0, 1 



i=0 



and the (exact) decomposition 



Gf^G 



(0) 



G'-p^ + G5J;" 



(139) 



holds. 

To proceed, we note that the summands of both S^, k = 0,1 are continuous and bounded functions on their relevant 
domains. Hence any singularity that arises does so because of the divergence (as functions) of the infinite series. In 
other words, the divergence is due to the 'large-^' contributions to Sk- To account for these, we can use large-£ 
approximations for the Legendre polynomials P£(cos7). We note first the trivial results that 
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where 

oo 

Tfc :=^(£+l/2)i/2-fee-*''P,(cos7), fc = 0, 1. (140) 



1=1 



Then we apply the Bonnet-Heine formula (see [18], p. 208): 

P^(cos7) = Pf(cos7) + 0{r^'^). (141) 

where 



P,(cos7) = y { ~ ^^^^^ ~ ^^^^ h + V2)7 - V4]| . (142) 

This relation holds for e < 7 < tt — e, < e < 7r/2. We note that, as above, the remainder contributes an overall 
continuous term to T/j, fc = 0, 1. It follows that the singular contributions to arise solely from 

00 

fk := ^(^+ l/2)i/2-^-e-^^"=P,(cos7), fc = 0, 1. (143) 

1=1 

For convenience, we introduce the notation 

to mean that A — i? is a continuous function. This gives a convenient way of indicating the removal of different 
continuous contributions to the various quantities encountered. We note also that in each case, continuous functions 
arising from infinite sums are identified by repeating the M— test/uniform convergence argument used above. In 
practice, this simply amounts to removing terms of order from the relevant series. That is, to obtain the singular 
parts of Tfc, we expand the non-oscillatory factors in inverse powers of and discard all 0{£~'^) terms. These do not 
contribute to the discontinuous part of the retarded Green's function. 
Thus gathering relevant terms, we have 

To = - ' ^;!2: {e'(^/^-/^)^(7-^.)-e-(^/^-/^)^(-(7 + ^.))| 
8 V 27r Vsm 7 L J 

+ l\f^^ |e^(^/2--/^)S(7 - r;,) + e-''(^/2'-/4)6(-(7 + ry.))| , (144) 
2 y TT sm 7 L J 

Ti = :^^/^^{e^(^/^--/^U(7-^.) + e-('^/2"-/4)^(-(7 + r7e))}, (145) 
2 \ TT sm 7 J 

where 



A{z) ■.= Y.^-\ 



B{z) :=^e^' 



(146) 



Recall that in order to obtain the retarded Green's function, we require the limit e — > 0+ of the terms above. This 
means that we need only evaluate A and B for real arguments: these have well-defined distributional forms. We give 
the results here, leaving the relevant derivations to Appendix B. For s S M, we find 



where 



A{s)^-V{s)+iU{s), 



4fc27r2 



V{s) = -21n|s| - 2^ In 

fc=0 

U{s) = -(7r-s)+7ri^i7(s-2fc7r)-^i/(-s-2fc7r) 

lfe=l A:=0 



(147) 

(148) 
(149) 
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and 



B{s) = --V'{s)+U'{s), 



with 



V\s) ^ 



-2^py 



2fc7r 



U'{s) = 



TT^5{s-2kTT). 



(150) 



(151) 



(152) 



Derivatives here are distributional, as indeed are the expressions themselves. For example, 2? is defined almost 
everywhere on M and is locally integrable, and so yields a distribution. 

Recall that our aim is to determine Gr, which is obtained by taking the real part of Gf, the calculation of which 
involves taking the limit e — >■ 0+ of A{z) and B{z). These functions are analytic on the lower half plane, which is 
precisely the situation that holds here: the arguments of B that we encounter have the form 

z — ±7 — rjt ^ ±7 — T] + i — , e — > . 

V 

This means that to determine G/?, we can simply set e = in the expressions above that combine to give Gp, and then 
take the real part (with the appropriate constant factor). The resulting terms are distributions rather than analytic 
functions. 

Of course we are not calculating the full Gr, only its discontinuous part which we denote G^"^. This is constructed 
from terms of the form ^(7 — 77) and B{'-f — i]) (respectively ^(—(7 + 77)) and B{—{'y + rj)), scaled by fa ctor s tha t 
include terms e*^^"''^/^ (respectively e''-'*'^''''/^). This can be verified by tracking through equations (137) to (1451. 



Determining the final result for the singular part of Gr is simplified by the following observation: the singularities of 
^(7 — 77) are concentrated at points with 7 — ?7 = 2fc7r, fc S Z. At these points, sin((7 — 77) /2) — 0. This factor mollifies 
the singular behaviour, and renders the corresponding term continuous. The possible combinations are essentially 
functions of the form sS{s), sPV{-), sins and sH{s), each one of which is continuous. The same holds for the 
singularities of B, and when the argument is —(7 + 77). A consequence of this, which simplifies the remainder of the 
calculation, is that there is no contribution to G^"* from terms with coefficients of the form sin((7 ± 77)72). Likewise, 



7±77 



E 

k 



A(7 ± 77 - 2/c7r) = ^(-1)''A(7 ± 77 - 2/c7r), 



where A(s) is any one of (5(s), PV^(i), In |s| or H{s) with singularities at s = 2fc7r. 

With these observations in hand, we can determine G™^^ (for convenience, we do not include the factor i7+(x,x'): 



R 



-l/2_ 



1 



7-7/ 



2kn) 



7 + 77 



PV 



1 



1 



7 + 7/ + 2kTT J 8 



cot 7 - 



cot 7 




V) 



In \ri + 7I 



00 

E 

k=l 



In 



1 



(7 + 77)2 



4fc27r2 



(153) 



The following observations yield the final result. First, we incorporate the comment above regarding the factors 
cos((7±7y)/2). Next, we apply a similar observation to the term 7/"^/^ go that 7; — !• 7+2A:7r in the 5 and U contributions, 
and 77 — >■ 2/c7r — 7 in the PV and In contributions. Note that in the series ^ 5, there is no contribution from terms 
with fc > 1, for then the argument is always negative. We then apply 5{s) = 5{—s) and retain the appropriate sum. 
A similar argument leads to a truncation of the series implicit in U and the X] term. Factors of the argument of 
the logarithm can also be removed. The coefficients of the singular terms in the series can then be identified: see (21 ) 
and ( 23 1 . Carrying out these steps yields 



fc=0 



t/(7 + 2A:7r) 
7 + 2fc7r 



(5(77 - (7 + 2fc7r)) + 7^0(7 + 2kTT)H{T] - (7 + 2fc7r)) 



E 

fe=i 



(-1)* 



U{2kTT - 7) 

2kT: — 7 



PV 



77 — (2fc7r — 7) 



+ zyo(2fc7r — 7) In |77 — (2fc7r — 7)! 



(154) 
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Next, we note the following results that derive from elementary properties of the distributions involved, and on sign 
properties of the arguments: 



1 



7 + 2kTv 



(5(r;- (7 + 2fc7r)) = Si-^rj^ + + 2k7rf), 



1 



-PV 



1 



2A;7r — 7 \ ~ (2/c7r — 7) 



PV 



1 



ln|r,-(2/c^-7)| = \n\^rf - ^i2k7: - 



Apply these in (154) yields our final result: 



(-1)'= 



fe=0 



fe=i 



t/(2fc7r-7)Fy 



1.0(7 + 2fc^)F(-ar")] 
1 



3dd 



i.o(2fc7r-7)lnK<5| 



where 



1 7— >-7+2A;7r i 



3dd 



(155) 

(156) 
(157) 



The designations even/odd are used to indicate that — q along a null geodesic that has passed through 



an even (odd) number of caustic points. 

We emphasize that the result (155) is exact in the sense that 



Gi 



where G^f™ is a continuous function. The result applies for values of 7 bounded away from and vr (see the comment 
immediately following (142 1. As noted earlier, caustics form at spacetime points with = 7 + 2kn,k > and 
r] = 2kTT — 7, fc > 1. The result (155) thus reproduces exactly the four-fold singularity structure discussed in the 
introduction, yielding both the four-fold sequence for the sharp term 

S{<t) PV (^^^ -S{(t) -PV 

and the tail term: 



Sia) 



e{-(j) ^ - In \a\ -9{-a) In |ct| ^ 9{-a) ^ • • • 

Furthermore, this result yields the coefficients of the singular terms throughout the spacetime. 

Figjs] shows a plot of obtained using the ^-mode sum Eq.(133). The four- fold singularity structure 
(5((t), P.V. (I/ttct) , — (5(cr), — P.V. (l/vrcr) , (5((t), . . . is clearly manifest in it. The £-mode sum is also compared with 
the Hadamard form in the normal neighbourhood. 



VI. DISCUSSION 



In j5], the following heuristic explanation for the singularity structure of the 'direct' term was given involving the 
Feynman propagator. The Hadamard form, valid for x' € J^{x), for the Feynman propagator is [2] 



Gf{x,x')— lim — 

e-i.O+ 2lT 



U{x,x') 



- Vix, x') In (o- + ie) + W{x, x) 



(158) 
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Gret - 
\ ^ 

(5(o-> 



-3 -2-1 12 3 




FIG. 3: (Colour online.) (a) (Left panel.) Visualization of PH spacetime. Here we have y = y' , the time coordinate T runs 
along the vertical axis and the angular coordinate 7 runs along the horizontal axis (which also includes —7). The dashed black 
vertical line corresponds to 7 = tt/2: this is a static path. Black dots on the left and right boundaries correspond to caustic 
points (for which 7 = 0,7r), red dots along the dashed line to points where null geodesies intersect the static path at 7 = n/2 
and the blue dot on the 7-axis to the base point x. Diagonal lines correspond to null geodesies; these alternate in colour every 
time the singularity of Gret changes its character, which is indicated alongside with a ±(5(cr), The brown-shaded region on 
the bottom corners corresponds to the acausal part of the maximal normal neighbourhood, where Gr = 0. The green-shaded 
region on the bottom diamond corresponds to the causal part of the maximal normal neighbourhood, where Gr — V (the 
bottom two diagonal lines, which emanate from x, axe also contained in the causal part of the maximal normal neighbourhood, 
but it is instead Gr — US (a) along these lines). 

(b) (Right panel.) Gr{x, x') as a function of At for y — y' {= 0), 6 = 0, (j>' ~ 0, 9' — ■k/2. The blue continuous curve corresponds 
to the e xact G r obtained using the 'smoothed sum method' of Sec. VII. D [5] (such 'smoothing' is justified in [19]) for the mode 
sum Eq.| |l33| l with m = 0, ^ = 1/8, Imax = 800, Icut = 200. The dotted gr een curve is V{x,x') calculated using a WKB Taylor 
series as in 20 . The dashed black curve is uo -|- vio from Eqs. (23 1 and ( |25[ ). The time for the red dots in (a) correspond to 
the singularity times of Gr in (b). The four-fold singularity structure (5(cr),P.V. (I/ttct) , —5(a), — P.V. (l/vrcr) ,5{(t), . . . can be 
clearly seen. 



where the biscalars U{x,x') — A^^^{x,x'), V{x,x') are the same as those appearing in the Hadamard form ([l]) for 
the retarded Green function and W{x,x') - like U{x,x') and V[x,x') - is a regular and real-valued bitensor defined 
in a normal neighbourhood of x. Using Gr{x,x') = 2H+{x, x')Re {Gf{x, x')) together with lim£_^Q+ l/(tT -I- ie) — 
P.V. (^) — iTr5{a) and lime_j.o+ ln(f7 + ie) = ln|cr| + in9{—a), it follows that the Hadamard form for the retarded 
Green's is as given in ([T]) above, and is valid for x' e M{x). 

In [S] it was argued, using the transport equation along the appropriate geodesic joining x and x' which the van 
Vleck determinant satisfies, that this biscalar picks up a phase '— tt' every time that the geodesic passes through a 
point such that 7 = or tt. In the case of null geodesies, these are caustic points. The previous argument was given for 
a static, spherically-symmetric spacetime and this is also seen in [16 using a completely different method. Applying 
then Gii{x,x') = 2H+{x,x')Ke{GFix,x')) to the 'direct' term A^/^/(cr + ie) of the Hadamard form, the four-fold 
singularity structure |A|^/^ (5(ct), |A|^/^ P.V. (I/ttct) |A|^/^ (5((7), - |A|^/^ P.V. (l/TTfr) , |A|^/^ (5((7), .. . for Gr. then 
follows. 
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One could apply a similar heuristic argument to the 'tail' term yin(iT + ie): considering that V{x,x') = 
A^/'^{x, x')V{x, x') with V always real-valued (therefore V would pick up the same phase as U after a crossing 
of a point where 7 = or tt) it would follow that the 'tail' term in the retarded Green's function should alternate as: 

\A\^^^ Vei-a), - \A\^/^ V\n\a\ /n, - \A\^/^ V0{-a), | A|'/' F In |a| /tt, etc. 

Using the Kirchhoff integral method we have proven that the suggested alternation after the first caustic-crossing, 
both for the 'direct' and the 'tail' parts, is indeed what occurs. The large-£ asymptotics of SecfV] give a completely 
independent verification of this result and extend it to an arbitrary number of caustic-crossings: the calculations 
required to apply the Kirchhoff method to evaluate Gu beyond the second caustic become considerably more involved. 
The large-^ asymptotics, however, raise a question regarding the causal character of the retarded Green function, which 
the Kirchhoff integral method successfully resolves. The extension of the meaning of the world function cr, as well as 
that of the biscalars U and V, outside the maximal normal neighbourhood is cleanly given by Eqs.([|), (H, and psf . 

Both four-fold cycles - for the 'direct' and the 'tail' parts - together provide a full-account of the singularity structure 
of the retarded Green's function beyond the normal neighbourhood. They are proven here for the PH spacetime but we 
speculate that this is true for a general set of spacetimes which includes Schwarzschild. Both methods, the Kirchhoff 
integral method and the large-€ asymptotics method, are in principle applicable to Schwarzschild spacetime. Although 
the Kirchhoff integral method seems to be of difficult application to Schwarzschild, we hope to be able to apply the 
large-^ asymptotics method to Schwarzschild. We note that the large-£ asymptotics for the quasinormal mode sum 
used in [J should, in principle, be compounded with a similar analysis of the branch cut of the Green's function - see, 
e.g., |21j : using in Schwarzshild the large-i* asymptotics similarly to here in PH would instead entail investigating a 
two-dimensional PDE in terms of the time and radial coordinates. 

Finally, we note that the four-fold singularity structure does not apply to spacetime points with angular separation 
equal to tt. As shown in Section IV-F, the singularity structure is two-fold: 6 — )■ —5 d . This pattern has 

recently been observed numerically in Schwarzschild spacetime [22], as has a corresponding pattern when 7 = 0, but 
with PV{l/a) — >■ —PV{1/(t) —>•••. We have not considered this latter case, but it should be possible to do so using 
the Kirchhoff integral method. The presence of points with a two-fold singularity structure is also predicted in 
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We have 



Let (/) e Cf^'iR). Then 



Appendix A: Distributional Limits 



1. Derivation of (77|. 



<A['\cI>> = V2 [ 

JR 



-2e,2e] ^ 

V2 -^E^] mdz 



J2. 



2e \z zVz + 2e^ 
z) 



2e Zy/z + 2e 



dz 



=: < a1;\ </> > + <A^^l<t>> + < A^l4> > ■ 

with the obvious term-by-term definitions of A^^-* , i — a,b,c. By definition, we see that lime_>.o+ A^'J (z) — V^PV 
Substituting z —: 2ey, we see that 
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Then 



It follows that 



Likewise, 



so that 



Thus 



We conclude that 



which is ( 77 ) 



lim < A^;\,^> 

e— S-0+ 



€-)-0+ 



lim < > 

e->0+ 



1 + 1 

= 2\/21n(l + %/2)0(O). 



lim a['j^ = 2%/21n(l + %/2)5. 



/c 



0(2ey) 



-%/20(O) 
-2\/21n(l + \/2)</>(0) 



dy 



lim A = -2\/21n(l + V2)(5. 



lim A('^(z) = V2PV , 



Define 



2. Derivation of ( 79 1 



A^")(z) := ln(y^T2^ + V2e)H{z + 2e), A2(z) := ln(yi)iJ(z). 



We aim to prove the distributional limit 



lim A^,'' = A2, 



which is (|79j). So let G C5"(M). Then 

/•OO 

< A^'^ - A2,(/i > = / \n{y/z + 2e + V2e)(l){z)dz- I \n{y^.)(f>{z)dz 



-2e Jo 

° ln(^/^T2^+V2i)0(z)dz+ In ( ^^^^ + ^ 



(j){z)dz. 



With the substitution z =: 2e?/, the first integral here evaluates to 



^0 1-0 
2e I ^\n{l + ^lT^)cf,{2ey)dy + eln{2e) J ^4>{2ey)dy, 



which, for any test function (/>, vanishes in the limit e — !■ 0+. Note that 



lim / ln(I + ./l + y)^{2ey)dy = -0(0). 



(Al) 
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Making the same substitution in the second integral on the right hand side of ( Al ) yields 

+ 2e + \/2e 



In 



0(z)dz = In 



Since G C5"(M), there are positive constants and such that 

0(z) = for all z > 2C0 

and 

\(j){z)\ < for all z e M. 

Hence 



It follows that 



Hence for any test function 0, 



^lhn.^°°ln(i±^)0(2e,)d, = O. 



lim < A^'^-A2,0>=O, 

e-i.O+ 



which proves ( 79 1 



(A2) 



Appendix B: Distributional forms of some infinite series 
1. Derivation of ( |147[ ). 

By standard Fourier transform methods, we have 

00 ^ 

^risin(ls)) = -(7r-,s), se(0,2^). 



Thus the series converges on M to U{s), the 27r-periodic continuation to the real line of 

0, ifs==0,7r; 



u{s) 



\{tt - s), if < s < 27r. 



This function can be written as 



Similarly, 



^ (00 00 ^ 

lfc=l fe=0 J 



Er^cos(^.) = iln(^^^-L_), .e(0,2.). 
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See for example Eq. 1.441(1) of [T7j. Periodicity of both sides indicates that the series converges almost everywhere 
on M to the locally integrable function of the right hand side defined on M \ 27rZ. The series is therefore a distribution 
which we can identify with |2?(s) where 

V{s) ^ In 



2(1 - coss) 
= -21n|2sin(s/2)| 

= -2(ln|s|+f;|l 



4fc27r2 ' / • 

\ k=l / 

This arises from elementary properties of logarithms and trigonometric functions and the infinite product represen- 
tation 



oo 

sin X = X 



Combining these results yields (1471. 



2. Derivation of ( 150 1 



Formal differentiation of both sides of (147) yields (150): the key point is that differentiating the right hand side of 
( 147 1 involves the calculation of the derivative of a distribution which always yields another distribution. Thus ( 150 ) 



is a valid distributional identity. The results (151) and (152) follow by applying the distributional derivatives 



--{H{s)}^5{s). 
as 



Some rearrangements of the series are also required. 
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